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AG moted carlice, both a foutth=-ordet finite difference 
acheme ant a secomd-atdert Finite clogent tethoed ate applied to 
the system fot numerical aclution with a atid intetwal of t4*9! 
Latitude, Additional teeta with different geid intervals and 
boundary positions have beer tade and results will be 
desctibed. Under discretization, » becoms a colum vector of 
dimension 105 with the exception of n=0 and 1, because for these 
tw wavenumbers, there are non-sofn boundary values to he 
determined for one of the dependent variahics and thus 7+ has two 
more detrees of freedoe. Correspondinaly, A becomes a matrix of 
9 blocks, each having dimension 16235 for all m other than 9 and 


hk. The details of these nurerical methods may be found in 


Appendix c. 
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ALE ow ER see OWE Rap ome ee 2S wl tt OMe Conde at loa!) mote) & 
ARLEN af ibe Bgaom tte woh at iom fs tas. CRY, we hehe Tt nt 4 
f&,JeASoletiume of Whe By SECA Mocha OD a hater? of watlal les, 
They fepend, of Couette, CH LeACapal dy oe tue Apmhatilcs of Cte meade fy 
Buk ahas oh Rhee bu tSat yg domi Plone, te weet acal Sepatatucan whi 
VRS Nahe ECe, RR wane tate t ate Ee comput atiotial acheme , tot 
ERe [lattes we Gelect the Gh: of Be Gf PLitegetmec and Teragte clement 
echohee fot Comparison, 

Pefoge BECUSEEM™ Bho gemefal qgirvtal made} we may Rortiaj det 
Eke wiFe¢cte of changing the geoncete;, Team paje:tospaje to a 
Charme. , wad aldo G2 abl iat the effect of hanging tt urcat ioe, 
beoO en Watping Ly. Toe RhGe PUShRae we hawe MAL Cjed out genetal 
More | £9 fepfcacrtt 2a chat? Freon 80°C to GOS amd toa jnej) ude cot 
Cech cutvature, “omejdet tow the shanme} tmede)] form uf jlizitsa 
a si fple decom? og fee Finite Api feger ce acheme , Ln ceanplje of 
the affect of chanting Ye from 4* Jat, te PO* Jat. may whe seer fy 
Fig. tt. We Show wm thie Figure the profile it Jatitydte of the 
eigenvector Foe both an ewvtermal (srivalest depth of ©,% bry and 
am internal fequivaion! depth of .42 Ent tode for wave suymt« t 
| one. Poth vectors are xf the Mweshy type with few zeros. Note 
Phat there if a Vistinsuiehahle cffeet wr the evternal mode by 
| changing the teurc22ianw, an effect which ja sant evident for the 

imternal mode. This indication of ar effect, altoudh act 

| conclusive, suideetea that care mist to taken in choneing a 


teuncation interval and that suck = cinice may have ac impact on 


the charact-rictic struct rres of the aode] even when thosn 
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LARCHER L veghin al pov, pce cf Fhe mode. , tionte How Sat tly theac 
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obeervabioan that the effect 3f chanrdgints Nowtdaticg and /ot 
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vlaneticy waves. Thue we ca: gay that hoetdaty effects show up 
74 SHuntary nwes for fhe chan; toda) with as counterparte if 
kre gnie-fo-onie case, ard tat curwature effects hawe a 
pramoaunced impact an noes of the torsest planetary waves. Rased 
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an *iminial valve*™ grusiem provided Chat the parameter ein the 
Seentency*® Rakfis amd 2 set of ieitia data (yu, Yo) are given, 
Wot cavh ia dievifuatl newte determined ai a2 Aumerical solution to 
CLP) we DNeett the eigenfrequcicy O.) and the equatorial values 
oF Se ye ate Yee) yensi fuel yres ta YRatsabize an integration 
Fowatt he FER pode, $f, at the pode, the boundary values of 
Ee Aahe ical cl pen ake wRufed ate Rot Met, «© 18 adyusted by 
VewRon*e fect heed att the jtteg¢atian repeated, More details on 
ENL¢ Ofodess wan be Fagmd In Append)s 9. The results from this 
PHO tafe ate cralwated ta cétabliah the nature of each mode of 
he 2S Fegonece oust fone, 

Serntal Criteria tay be utiliee ta compere the results of 
Bho Tehinck pmg*® eteuctures with Ehoae of the fiaito-difforence 


etyatione. We yee the foallawing: 


ay fatia of Trequencics 
4) Loaet aguate Bifference of structures 


€) goetelation coefficient HWetween the structures. 


“he tatio of fretpencjice nay be cetabhlished ones the shooting 
croecetute Xae converge, This occurs when an adéustment to the 
frormrcy ferivet? from the difference syaten no lonaer changes 
she boundary calm at the mnie to wtich the shontina procedure 
insogatos, asd tmpecver, thre hountary vaiue corresponds to the 
heurdare vale agclicahle to the ficite-differenee solutions. 


ALPoerateiyv, whe the Mumdary valucs of the two solutions 


correspand £7 a oreeelecte4 «enall 4a fference, the shootina method 


iS ait 8 Rave Guofvet yet to a The we*™ pps scai made. df, 01 
Course, Che beutiaf x vatse .5 fel ab @uached ty Ube ahautiag 
hethost Before Ofte Bfenpeetcy chatiges to Om Lh the eiegihygty of a 
fiffefent cyjefiv¥alwuc, then wm Ray Slate that the Ep Alter 
Qiiterenhee BLEGELOFS Gh Gaeskiah 4s ihdeed a catpetatianal made ; 
bee, ER there MEE Conte spams Lo ang Made whsch satsafies the 
Hitfecential equations, The gatia al Tre quctc.cs thud measures 
the Converyed Creyuency vaso of the shaat .hg tethad ta the 
Atiginal Crequoney with which the methad tegah., Vajuca of this 
ratio tac Etom unity Cleatl, tuggcat that the finjtesd) fleretce 
ayater Made cither dooce sat apprasimate the tt ac tode weld, at 
that it Fepeegenta a Gonrpulatianal tode. 

tn copions of the salucionm whee frequctcy separation jis 
Large, 1.%., whete sequcntics Frequencies which gatiafy the 
equations change their valucs gubstantially, tte tatio of 
frequencies ney fot be a Gatiafactory neaaste ta identify 4 
cComputatiroral mode. Nder auch cotditiane if if NeccsPary to 
have aditional measures and we utilize the aonoari pon of 
Pidgenvector steuctures, Such corparieannc are nade coteurrent ly 
aith the catio-of-frequencics, inplyinag that tre shooting method, 
if not converged, has at least heen compicted, The two 
;rocedures identified above as 6) ard fol ativnive a comparisor 
oF the Latitudinal structircs at all poicts available from the 
chose. truncation. Since wr age %8 = 5° latitude, the 
comparisons are dqenerally tade for each physical variable at 3% 


MOints. Ceilizing these peints for both vectors [the original 


solution and the “shot” ssistion) we caiculate both the least 


Spwatese Mit ferceice awh the Coteetatin Cuethiclent, "ese 

Re sSuhee five whit iomhal LAs l yh 24h flaw Ete two aoletions 
CIMPate. f§ Ray Sapece hat thie CHeuducholies Combate wed! fat the 
SL sky kce RI ows, Rtetmatceby, the ee Wachcies may te dittete ht 
fin ERo@e Chetial fesplome alluded fo atvwored) while the steuctutes 
Ray Me BLNiLL ae. Witlalty, gf Whe Slt uct sted thawe sjgnlitjicant 

amp litQyphe im cAly a Shall Latitadjimal tatge suct at the 
eteahoe ial fegion, OQ2f Lali stical meangtes om Sttuctyte tay 
G@eill peove Ro be jmadespuate ho gice 4 def jnjptive andwet te the 
Reyec MaAkuee of Khe ALF fe gemee segya' ,or was yt jor, 

To domMangitatc § he APH i itatian of the aghast, gq tmethad and 
the patam@etete of conhartigom, we present Tabje @. $8 this table 
we Compare thn def af weatwatd prapatating grarit waves whose 
SEfuckutee feakablithed by wohetion of 1]70) ate weer phe d on 
Tige, FO amt Tl. Me deceit] anlytiong For the bt order 
*itferonce eqyationé ioe plavctary wave fie ate’ fog the catetnat 
verfical tewie., Phyaica: as well ae cotputatianal modes are 
ime luce? pm ERS amt feecleding Sesty todee and «© actward 
OTIWATIL ims teacity fevde gs aud jpadeet without add t jaca) 
eriftetia, we would come der all of theaqn modes ac physical. tet 
ue aw determine what thm shanting method atrlics aboyt this 
aot. RA the domain of tow inde ew vabose  ige., 22-381 we ger lope 
values of leaet ertaates couried with nedlidible valuyce of the 
caerelation cmefficiernt. Trees mwtec are untnuttoediy 
comoukt ational amt their etenctsres fia, 194 surnoort this 


come lueian., Ngte that the ratin-nf-frenueccics pt at a 


coneiagtestive dort indicator for calecticd cormnytatincal andec, 


f3nte #. DEabESEiSS oh Feo al tS with Ule af aat. ng he T Had (qs. 
L¢) appbiec? to Case for a-8, 659,595 Am, @tn afdey differencing. 


Mere > fate f bh OGsle G of Ugedqueehecy value 
taba = Aafbe@ OF sof O.<CFOSaiNYS TA StEuctue (4: taa high) 
c= Fecogehiey SObve d EG an made. caveat ivns 131 ~3074) 
ra - feyqucicy Rad. fod ty shast.nag method 
mete = leaQt aguefe Apffcrenece af nade} and ahat Bt ucture & 
Sug et vorgelation cocfiteieat af mode} and shot structurcs 
So, trdes iz) fn fasts + c.c. 
a) * 7.99 #109, =13.7 41.6 -0,07 
ia tl ?,e8 8.0% 7.02 0.52 =0,4) 
i * 7.69 #04 1.0% 0.40 =0.09 
4 0 7.49 *.44 06.98 0.41 “0.21 
\e * 2.26 7.44 t.Q02 a. 5 =0.00 
Qn 9 7.24 4.%6 0.79 0.33 =0.23 
2 ” 6,87 §.23 0.76 o.23 =-0.16 
e2 * 6.74% & 86 Y.02 G.31 =0.608 
2} ? &.46 6.86 1.06 0.91 0.95 
24 * 4.16 6.30 1.82 G.279 «0.00 
“% 6 6.0} &.30 1.64 0.Aj 0.986 
a6 5 $.59 5.76 1.03 0.59 0.99 
a” ~ 4.5) $.76 }.05 0.25 =0.10 
a? 4 6.43 $.2) 4-02 0.009 1.90 
«? * 4.78 £.23 1.09 G.25 =-0.00 
to 3 4.67 4.73 1.41 a.00 1.46 
VW 2 4.20 4.23 1.01 6.90 1.60 
Ve . 4.00 16.8 4.21 0.20 -0.12 
rr ! Le | 3.74 }.00 0.60 1.00 
1g oO 3.24 3.25 1.90 0.80 1.00 
8 + 3.13 3.25 1.04 0.05 -0.01 


For mate 624 the ratio is only two percent tram unity whereas tor 


632 be t6 a9 factor of four. Correspendingly, acqghiqitly ema)! 
vatues of least suuare ditferenees ate coupled with unit values 
of tne correlation coefficient, inficatiog true physical modes, 
for these nades the catio-of<frequency 38 also reasonably clase 
Rear UALlyY, Gubstantiating our interpretation of theae physical 
Aatu@ce, 

The shooting method has been applied using a 5° increment 
for integration fram equator to pole, a value which corresponds 
tor the increment used to solve the original dificrence 
ritustions., Recause we felt that such an increment was too 
cnarss, we tested the process with much smaller intervals, the 
smallest representing 0.01 deqrees of latitude. The results of 
that experiment were striking. As the increment was decreascd 
from five degrees, the solutions deteriorated, showing the worst 
cesults for an increment of one dearce. As the increment was 
fucther decreased, the results improved qradually. We concluded 
that the five degree increment for which we gave a demonstration 
was satisfactory and could not be sicnificantly improved upon. 

For the high index modes (Nos. 15-21) our statistics do not 
indicate any physical modes. However, Fiq.10 suqaests that sone 
of these modes may by physical. To climinate this ambiauity, we 
have modified the shooting method, hopefully to make it more 
sensitive to high index solutions. Cather than reducing (17} to 
two first order equations as was done in (18), we reduce the 
entire set to a single second order cquation in one dependent 


variable, say + . This equation may be expressed as follows: 
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(19) 


,@ ; 2 ee 
2,2 *. mf zt 
Cle.) & 2 (wv =f") = wha ep oees C) -e hos 
c? cos’; reine vente 


Urilizing (19) as the shooting cyuation, two initial values are 
requiced. These were chosen as the values of » at the equator 
and at the first point (5° iat). For this modified procedury 
various increments were also tested, but again such variation had 
Little impact on the solution. 

Results for high index modes using this procedure are 
described in Table 5 which 18s identical to Table 4 except for the 
shooting method applied. It is evident by noting the correlation 
coefficient for modes 16, 1a, 20 and 21 that the new technique is 
considerably more sensitive in this range, and that the 
correlation coefficient, as well as the other parameters, 
identifies these modes as physical rather than computational. 
This effectiveness of system (19) has been tested on other 
planetary waves and vertical modes with equivalent success, as 
»@ll as for the finite-element method, the results of which are 
described on Table 6. 

Consequently, we have ietermined that the appropriate 


shooting method to use is based on system (19). We have checked 


caretully to establish that this procedure also successfully 
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Table 5. Statistics on results with the modified shooting method 
Fq. (19). All other conditions are the same as for Tabcle 4. 


wr eee ee = = wineedee Be Sete, 2 7 Se, 4.4.86 = or eee 1 


No. Index ts fm (m/fs eis CLC, 
ll ® 8.47 8.37 0.99 0.57 -0.49 
12 bd 3.46 9.47 .J2 0.5) -0.33 
13 by 8.17 86.06 0.99 0.85 -0,47 
14 * 8.11 7.64 0.94 0.69 -0.29 
15 bd 7.99 9.12 1.14 0.43 0.23 
16 lt 7.88 8.06 1.02 0.22 0.95 
1? ° 7.69 8.79 1.24 0.75 -0.08 
18 10 7.59 7.66 1.0) 0.11 0.99 
19 bd 7.27 7.49 1.03 0.90 -0.03 
20 9 7.25 7.27 1.00 0.06 1.00 
21 8 6.87 6.85 1.00 v.03 1.00 
22 * 6.75 3.72 0.55 0.87 -0.02 
23 9 6.46 6.43 1.00 0.01 1.90 
24 ad 6.16 8.27 1.34 0.87 -0.06 
25 6 


4.03 6.00 0.99 0.91 1.00 


me ae ee ee 


Table 
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Same a8 Table 


“xcept tor finite-clement expansion. 
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0.9) 
0.57 
-0.01 
0.01 
-0.02 
-0.41 
0.32 
0.99 
-0.06 
0.99 
0.99 
0.99 
-0.04 
1.00 
1.00 
-0.06 
1.00 
1.00 
1.00 
-0.01 
1.00 
1.00 
1.00 
-0.00 
1.00 
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idencviffied the computatioral modes ir the lower index range as 
was demonstrated for system (18). Te definicively identify a)] 
possible computationa]) modes we use all three of the parameters 
specified: (a), (9), and (c). Where possible, we also use the 
index parameter since when the zero crossings become too 
frequent (~24¢) the solutions are clearly computational. 
However, such a test is best done visually, and as such iss not 
only subject to human crror, but is ¢rohibitively time 
consuming. Thus we use the correlation cocfficient as the prime 
Parameter, noting that unless its vaiuce is greater than .95, 3)) 
indications are that the representative structure is 
computational. If the condition is met (c.c.>.95), the other two 
criteria are checked for consistency. In this way we can 
definitively identify and isolate al) computational modes. 
Finally, there are a set of modes which cannot cxist 
physically and must be computational; these include all the 
eastward propagating Rossby modes. They may be identified 


immediately from their frequency and the shooting method 


corroborates this identification. 


2 


5. Physicel vi. Computational sales 

Having established 2 Sotho whereby the distinction tetween 
physical ant computations! stoden characteristic of Gur Model is 
Mate evident, We may ask about theie relative distributions. 

This ssuestion docs net appear to have a utique anawer, but BOM 
tntormation if available (roe our ca.culations, a8 indicated by 
Tables 4, 5, and 6. We note that tor the external vertical) mode, 
the longest planvtary waves have the most physica! anodes. This 
relationship 18 true alse for the internal modes, hut there is 4 
Ceduction in the number of ohysical modes as ane procecds to 
higher internal modes. 

Consider the cxample of westward propagating gravity modes 
for the 4th-order difference model. For the external mode, waves 
zero and one have half phys.cal and half computational modes. As 
one proceeds to the shorter scales, the ratio of physical to 
computational modes drops asymptotically to about 20 percent, 80 
that by wave m=#1l3 this minimum ratio has heen reached. This 
minimum ratio seems to apply also for the first three internal 
modes although the ratio for the longest waves (m20,1) decreases 
rapidly, so that for internal mode three (k*4), wave number one 
has only a 25 percent ratio and drops to 20 percent by wave 
number nine. 

We note further a subs:antial reduction in physical modes 
between internal vertical modes four and five, and for yet higher 


internal modes we sec only three physical modes, even at the 


longest wavelengths. The results clearly indicate that for 


x 


highly structured vecturs 8st et maves ast Cate ral modes} there 
are Rany Nore Gath slatiotal Greate ihohienks in the fatal of 
‘Seavity waves. 

A @imitare felatiomehig ee Ble tee Ue hosat sy trade s alttiaug!t 
more pRhyeiea!l modes Fetal Tar the sturtet wave fenaths Ctar. 
appear for the gravity tertes., We do fet have sufficient Gata to 
Comment on the increage of Cumputaticfal khoashy mades witt, 
increasing interial mxice, byt the Catfe}ation docs cayst.,  ladeed 
it i@ difficult to tdentifiv any true physical modes assoc, ated 
with any wavenumbere for the Jaet (ninth) internal trade, 
Although there if seqliqible atplitude aagodiated with thie made 
in atmogphecric data, Such coapytatioral econtratutiane tay wed] 
infect and deteringate nantinese inteqratjions, 

te is of interest to mote that "he physical modes deve oped 
from the model equations, at least in the gravity domain, are 
comarkably similar for the 4@th-order and finite-clement 
representations. Perusal of Fias. I and Tl take this compar isin 
clear. The computational sodes do show some atructural 
‘differences. This can be seen hy corparina the enavelopen of 
comoutational mdes for corresponding vectors on the two 
fiqures. Tf all sodes are used in nonlinear integrations, 
because of the differences in the corputational modes, the two 
systems should yield different results. If, however, the 
computational todes are filtered from the data ard only the 
physical modes represent the initial conditions, ‘he two systens 


(4th-order finite difference and fin’ te element) should yield the 


Same intedrations until the cortmrtational mode atoclitudes deve lon 
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H, Emi kiali sation 

tm the previous ae Ct ios me hawe tr ghiighted the method 
whereby computati.o al tievtes tay be separated Trot physical modes 
when Characteristic Sr ati ats oF a litecatived syaten ate 
deteenined., tt has align beon demonstrated that if a nonlinear 
Mote, 28 expanded in theae characteristic modes, filtetitng of the 
imitial data may Save a sidsificant aapact om the integration 
properties of the system. im the past such filtering has beer 
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“hat vechaps computational crrors dee to rumerical proceases 
(Finite-differencind, ete.) nay aisn te ashihited by filterina 
initial data of unwanted comritatioral modes. 


For complete elisination of computational modes, their 


initial values as well 4s higher derivatives i- time must he 
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Tero ', tecreecnt the latitdinal arid noints and c, represents 
Rhe nine och? eigta feveie, “Se vectors G, are those which were 
terive? in appecdie Ro and deecrirmed cava Tid, 2. The amplitude 
coefficier*s uy, we to he detertined. By cormalizina and 
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The process 13 sore simply ccpresented ft, 2 matrix equation if we 


detine the variosbles a8 tol Own, tat te @ tat eix with 


op 
elements Unie? Of dimension (4749) and jet cach of the 

vertical structures be represented by a vectur af Aine 

elements (Gytol)) or Gy If we now develop a matrix G which is 


made up of all nine of the vectors G_, then we have 


x’ 
UF Une (23) 


where vn i8 a matrix cach of whose vectors represents the 
projection onto the corresponding vector of G. The latitudinal 
structures of these vectors correspom to the profiles which are 


analyzed by our model. We shall denote these vectors Ua where, 
Coe het (24). 


Tt should be evident that there are nine such vectocs for each 
Fourier coefficient m, and that there are in addition an equal 
number for both the v and » fields. By stacking these vectors in 
the order defined by (20) we: may writs a matrix of all variables 
us Ane Note that for each : and each & (vertical mode), the 
vectors of te are created b’ first setting down all latitudinal 


point values of Une followed by all iv, and finally by all «w/c. 


™m 


The total length of these vectors is clearly 111, and they may be 
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6.) faturidinal projection 

At thie 6tayjye of the transtorastion process, we have our 
4ata available for each planetary wave (m) and for each mode) 
vertical mate («) a6 a vector on the latitudinal arid, to 
eagtablish the contribution of this latitudinal profile ta each ot 
the model modes in latitude, we siaply project the data vectors 
(which we naw have) onto those eode! sodes. Rut those mode) 
moxtee depend on the numerical process used to represent the 
model, SO wo must froject the data on both the modes of the @th- 
otder system as well a8 on thoae of the finite clement system. 

These model vectors have been discussed in Section 1 and 
Appendix C. They ace the eigenvectorr of the matrix & (ace Fa, 


17) and nay be defined froa 
Azsa svt (25) 


where § is made up of the vectors 5 These vectors are 


mr” 
displayed graphically by Fig. 10 for =26 and kel (external mode) 
yased on 4th- order differencina and on Fig. 11 based on finite 
element analysis. Protection of our ‘ata vectors in terms of 5 


yields, 


(26), 


and finally, the projection astplitudes are determined by 
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inversion, Thus we can establish the relative strength of cach 


Of the model modes Sn , 19 vur data sample from, 
e 


fad 


oo 
* 
LD 


m,k Am k (27). 


It must pe recalled that $ depends on the numerical method 
applied and that we have two separate matrices to consider. 

We are now in a position to filter the given observational 
data of contributions from computational modes. Thus, if the 
filtered data is utilized in a nonlinear version of the model for 
which the filtered computational modes are representative, the 
numerical integration of those filtered initial conditions should 
show, hopefully, reduced computational effects. Since each of 
the elements of mak cepresents the strength of the corresponding 
vector in the data set, let us first establish how much amplitude 
of the data is involved in the computational modes. Recall that 
we have determined in Section 4 how computational modes are to be 
defined. 

To show how these amplitudes distribute amongst the modes, 
we present Fiq. 12 which is representative of wave m=6 and the 
external vertical mode, k=l. The corresponding structures of the 
AOdes are presented on Fig. 10. On Fiq. 12 we have connected the 
physical modes by line seqments. For this case of given (m,k), 
41 of tne 105 modes are physical and the remainder are 
computational. The amplitude squared (ee a) of the 


computational modes when summed accounts for only 8.8% of the 


total squared amplitude in the data. Since the physical 
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variables represented are flow and geopotential, we may also 
equate this squared amplitude to normilized energy. That the 
physical modes dominate is evident from Fiq. 12. This will 
hopefully be the case with the other scales (m,k). A similar 
distribution for the finite-element method may be seen from Fig. 
13. 

We now reconstruct the latitudinal profiles of the data hy 
inverting Eq. (27), however discarding the amplitudes of the 
computational modes. This is done by using only those values 
of em which represent physical modes in making the 
calculation. Clearly in the example cited, 8.8% of the squared 
amplitude of the data will be removed. These reconstructed 
profiles are now in a format to be used as initial conditions for 
numerical integration with our model. Fig. 14 describes how 
filtering alters the original profile. The figure describes the 
latitudinal profiles of Un, Kb) both before (dotted) and after 
(solid for finite difference, dashed for finite element) 
filtering of computational modes. The integrity of the profile 
is essentially maintained, but extremes are smoothed. MThis 
indicates the impact of the 4th- order finite difference operator 
as well as the finite-element operator on the profile and how the 
computational model is able to deal with extreme latitudinal 
fluctuations which are actually observed. By filtering out 
unmanageable fluctuations, perhaps these factors will not impose 
strongly on the nonlinear progression of the solution. The 


comparable results hetween the finite-element system and the 4-th 


order system exist because their physical modes are so similar. 


Re vben- _ wyaben es a - ager eS fo 
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7. Nonlinear integrations 

The computational modes identified in the previous sections 
have their impact during a nonlinear integration of the equations 
of motion. To assess this impact, we propose to integrate the 
nonlinear version of our linear system of equations as first 
Presented by Eqs. (1). Our intent is to define the effects of 
latitudinal space truncation, so we may represent our equations 
in terms of planetary wave numbers (spectral expansion) and in 
terms of the model's vertical modes. 

Since we have identified the vertical modes, upon which the 
wind and height fields depend, with the vectors Gy, (9), we may 


expand those variables in terms of the G vectors; 


<i 
t 


= : VS, (6) (28) 


y= pG (a) 
k k”k 


where the vectors G, satisfy the equation, 


oo (de 28k, aeck 

36 2 90 2 

and 

: Gy (a; ) Gyr (o;) = Sys (29). 


The prediction equations as nonlinear forms of (7) may be written 


as 


oh Pe 
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If we neglect vertical advection, (30) represents a set of three 
nonlinear equations to predict V and y. We may now substitute 
the expansion in vertical modes given by (28) into Eqs. (30). 
Utilizing the orthogonality of these modes, we will multiply by 


each mode in turn and sum over all points in the vertical. This 


yields, 
av, OD oe pie es = 
Be i ate sy Gan RG Rey er Ke 
v 
Ivy 7 eee y 
Be ea, pia Reiger ee et (31). 
v 
OK K', kt = : Gy (9; Gyr (9; Garo; ) 
2 2 
B yp cl ei Ven a ne anueny: < “k 4 
kk ket 2 © x SFG (OIG ar bos) F my Oy pr en 
1 Che Che 


It is clear from (31) that because of nonlinearity, the vertical 
modes are interactive and that the amplitude functions 

Vy and ¥, must be calculated for all k at any time before time 
extrapolation can continue to the next time level. To simplify 
the calculation but still maintain the impact of non- linearity 
on the evolution of the latitudinal computational modes, let us 
reduce the system to a "Shallow-water" one, by simply presuming 


that only one vertical mode exists in the expansion of (28). 


onder this constraint the summation on the right hand side of 


(31) vanishes and the coefticients become, 
a, = 8, = 2 Gi(c.) (32) 
k k i kt i ° 


Since we shall now consider each vertical mode {k) individually, 
no loss of generality will be suffered if we drop the "k" 
subscript. Let us now refer back to Eq. (16) where we have 
expanded in longitudinal waves denoted by the wave number m. 
Furthermore we can recall the vector of variables x as it is 
defined for each wave (m) and one vertical mode (k). Using the 
operator matrix, A, and noting the orthogonality of the Fourier 


functions, system (31) can be rewritten as follows; 


—_ u 

3 : _ = 

set iA(m KOK = Gg, = (33). 
4y 


This equation applies for each wave (m) and mode (k). The non- 


linear vector on the right hand side of (33) is defined as 


follows: 
ima u du voau uv 

Gy 2 a, fe ™ al sates ay + ¥ xy - EE tans] 
@. = <téy J -imd, ( u OV ee VOBY 2 u?tang) 
ly . k acoS@¢@ 9d a a9 a (34). 

TOK -imd u ap ov ay 
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It is an easy matter to remove the A-dependence from the q's by 


rr rere 
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integration. Because of the nonlinearity, as we expand V or y in 
m, the q's will generate a double sum over all wave numbers. 
However, orthognality requires the simple addition rule to 

apply. Therefore, the double sums are reduced to single ones 
over all wave numbers, since we have that m’'=m-m'. As an 


example, 


-a, fe7im u au dvi = -2ni £ (35). 


acos¢@ or ~ acosé “k m' 


Eq. (35) represents the first member of q,, as described by (34) 
and shows clearly the nonlinearity. It must be noted that the 
summation goes over all allowed values of m', here chosen 

as |m'|<20 and that m in (35) refers to the wave number (m) 
associated with the vector x in Eq. (33). 

Although it is now possible to expand (33) in terms of the 
latitudinal modal functions developed for the linear problem, we 
shall simply convert the system by numerical means to a 
difference system using (a) 4th-order differencing in latitude or 
(b) finite-element differencing. Both these procedures have been 
described in detail in Appendix C and need not be reviewed 
here. However, it should be evident that the left hand side of 
(35) becomes a matrix set over all latitudinal points and would 
be satisfied by the normal modes (including computational modes) 
if nonlinearity did not exist. Thus the vector 0 must be 
expanded at all Latitudinal grid points and the appropriate 


difference operator applied. Because of (34) and (35) we sce 


that all grid points and all waves interact to affect any wave at 


ee oe 
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a particular yrid-point. 

Given an initial state represented by x(t=0) for all grid 
points and planetary waves, as well as for the specified vertical 
mode, the value of xit) can be determined by a suitable time 
stepping procedure. To accomplish this, we utilized the well 
tested leapfrog scheme, using a multiple forward start. 

At this point the development of initialization as discussed 
in the previous section comes into play. System (33) in its 
final form for integration clearly involves computational modes; 
indeed we have isolated their properties with great care. We 
have also, however, removed the contributions to these 
computational modes from the initial data. Thus the integrations 
should begin without computational modes, but it is clear to see 
that because of the vector 0, these modes will he regenerated 
during the integration. How rapidly this happens and how 
strongly it depends on the numerical form of the equations will 
be established by integrating system (33) for both numerical 
methods with and without the filtering of the computational 
modes. The growth of the modes can clearly be identified if, 
during integration, they are monitored. Thus we will be able to 
establish their impact and the benefit of filtering 
(initialization). Finally, if growth of the computational modes 
is rapid, they may he Filtered periodically during the 
integration thereby inhibiting their effect. 

Computer programs for the numerical integration of (33) have 


been prepared and checked, but no successful integrations have 


yet been run. We cannot therefore report on these experiments, 


4s 


although we are actively workina on this project and hope to 


complete it soon. 


8. Conclusion 

In our search to isolate and define computationa) errors in 
atmospheric modeling with intent to cefine the “best* 
computational schemes, we have uncovered a method which, although 
not yielding the universally best nunerical scheme, will 
nevertheless provide the needed insight to choose the proper 
scheme for any particular model]. By model we mean here the 
complete set of model equations including all physics, dynamics, 
boundary and initial conditions. The procedure is oriented 
toward space truncation and we have focussed on truncation in the 
horizontal rather than the vertical. The essence of the 
technique is to define the normal modes of the linearized version 
of the model, to expand the initial conditions in terms of these 
modes, to separate the physical from computational contributions, 
and to filter the computational modes from the initial data. The 
model is then integrated with the filtered data and refiltered 
during integration, as needed. 

The primary issues associated with this process are the 
selection of normal modes for the linear model, the 
identification of computational modes, and the filtering of the 
initial data. We have described thes» procedures in detail by 
application to a specific model and by using a high quality 


global data set. The model used is tised on the current non- 


forced GLAS global model usiny o-conrtinates in the vertical. 


Th: mostel os litcesactzed on a S€ate of gest @fAt 3ts vertyca. 
dependence is Ccansfucmed C9 nOortal addes Ly Converting to a 
Seoucated adblteconee system asing second-order differences. 
Conventianal boundary cond). pons are appiiedt at the top and 
bortom of the model. These conditivonas could have gome (ai thous 
probably not too significant) impact on our ecunclusions, and more 
renlistic top boundary conditions are discussed in a companion 
report. The implications of vertical differencing as chosen here 
also deserve additional study, but again, variations in that 
procedure may not have a dramatic effect on aur findings. 

The horizontal equatiors, a8 separated from the vertical 
ones, have themselves been separated into longitudinal and 
latitudinal modes. Since we are here interested in global 
prediction, we have used the periodicity propertics of the 
earth's atmosphere in longitude to represent that coordinate by 
Fourier series, thereby separating out a set of equations in 
latitude and time for each planetary wave in longitude. This 
process inhibits computational errors due to differencing in 
lonqitude, but aliasing errors may still have an effect in a 
nonlinear calculation, due to series truncation. The normal 
modes in latitude may be established once those equations 
(depending as indicated on ooth a planetary wave number and a 
vertical mode) have heen converted by some computational 
approximation to numerical form. 

We have used two independent computational methods to 


transfora these equations; a fourth-order finite difference 


scheme and a finite-element method, both on equal jrids of five 


SS a —— 


degrees Latitudes feam the funth to the south gule. WM uss tw. 
AiR ferent methots we ate able Ga ime tify tte diftegences t'at 
fumerical 6ehenceé show whem used ta integrate vumpdea floes 

eCHuat tons, AS capected the fuettial medes (at develdoal team these 
CQual lOnS Separate ints pakert lar gpas ity and hossty mades, 
Information from 8*hie procedure has teem used widely to tytees 
LAER Lal data of unwantat Freueney Cofarsty! wave ptopadatins, 
THe Geparation of the gravity “Raseby Mades foe the Made) af sug 
CHOLCe Ghows youd comparability with atthe made} made sefatat oot 
Cigeueged in the Literature, Whe ape cific peaye of ug 
HiscusSion cefers to the dethinet jon jn this modal @et fietweer, 
computational and physical noades, 

Computational modee are Cloatly identifiable undet ceptayn 
CaONditions, almost impossible to jdentyty gnder others. One rust 
therefore egtahliah a systenatic procedure which will] identyty 
them. The technique whieh we have evglved teatea all moter 
developed fram the nurterical wodel te determnine af they satiety 
the differential ectuations al@o. tf they do, they are true 
physical modes; if not, they are computational. Jere cross, ose 
of modal structures also dive tood itsitht isto the dietinct .on 
between computational and physical ew des, but may not he 
conclusive. Based om three criteria identified and defined as 
our teport, aS well as chonsind the -y6t succeasful procedure for 
testing modes in the differential escations, wo believe wo car. 
separate the ohysical From coenrutaticnal modes of a numerical 


model, and do so definitively. 


To show how the co*putationral ~ des itnact so me} 


eho fh ah poh Bute aah eg) ow toate Uta Seb be Pte tect a gaalba. Gat. 


teh On all Rthe fete D Reade o Beet cope Tot oot Teet Made J ata 
Salcwbake fhe Softy aki otal Cate itul,sans, Te Vata sect was 
RaecS Ef oh the isn, sO), ahalyeged iy tioWar, Me Feau]ts ative 
HR anh Foe 38 be ast ute wate etl ree ecay thante, PLA Pref Ce ht of ttc 
Neg &f4eate abel obec of Rhee Bala Peotected ahtu the 
Compa kakponal MoteS, ALE ee This je tot 2 tegligjtje 
afp be budte, FEHPEGtCEE OM of Ihe plz Sical tode data (aftcs¢ 
EF, Lech ime Ehe Comput at poral modes) yicided a prafije Far the 
WALA Lal haha MOE Suhel antiaily altered Fram the sonfjijtereda 
PEO Le Bet WL EN coewhat fcRuged cet eencé, Wie Lanai ratce that 
fhe FehRefing Process Can FoMawe the jApart of AAnpulat jonal 
1I_~ee WLEROGE Recceearily Fotoring the phyei¢al ecantent of tte 
tata, \ 
The Final step pm tegtitg the ef facjency of our technique jn 
Pr arte state fhe meme ri cal oduatiomea with the Fpltered initial 
eoab pb bome. Thee procedure wil) tetoenite how fast compytational 
Rates temw Fe ENtograt pw and tow they inpact on the 
mom inear cvebutionm of the aniathiom. Moreowet, hecause we tave 
chosen twn separate como stat iota, actenes, the effect of 
Pranear yan Fype om Fhe tte teatior peoveess wii hecome evi tort, 
we Nave autliae? the orocedete for compictinag this last ster, and 
we have prepare} the proarais rweessaty for its implenentatior. 
nioetunately, theee calculaicss hawe cot yet been completed a4 
we cannot report am their ontecm. We are continuing this 


eeorgiment and will report *he res.its as soon as they are 


available. 
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fe tS our contention that the procedure identified and 
agutbiaed iA this report should prove valuable to the prediction 
community by reducing the etfeets of difterencinag 
approsimations, Mien our test integrations are completed, we 
hall apply the procedure to a complex, forced model which is 
Cugtently operational at CLAS. Finally, although our efforts 
have focusted on global sodele, we belicve that with some 


netifications, our Scheme can alfo t appiicd successfully to 


Limited area nedele, 


Appendix A. The thermodynamic equation 
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Linearizing but retaining the vertical advection term, the 


thermodynamic equation (5) becomes 


From the Poisson equation, 


Py R/Cp 
0 = T(p-) (P, = 1000 mb) 


and using the equation of state 


Pa = RT, 
one has 
R/C, 
Q = p°v/©p pe ers 
and thus 
aénO _ Cy l ¢ 1 aA 
“OS C_ oa A aa 
Pp 

gino | Sv Lan , 1 da 
“ot! Cc” WF ote A dt 


(A-]). 


(A-2) 


(A-3) 


in which perturbation quantities # and a are negligible when 


compared to quantities of tne basic state WN and A, 
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respectively. Also, from Fq. (6) and the perturbation equation 


of (4) one obtains 


2 
day 1. a4 3A On 
se "9" sexe * oss oe? (A-4). 


Finally, by substituting Eqs. (A-2), (A-3) and (A-4) into (A-1), 


one gets a thermodynamic equation with the following form, 


or 
2 
a ¥ 2 «Nn? - - 
5t30 NU w 2 -Bs (A-5S) 
where 
Cc 
203 ee vii 1 aA 
Nv” = B vee . + kK 35) 


Appendix B. Solution of the vertical problem 


The second order finite-difference equations for solving H 


are 


1 
——z (-2H, + Hy) 


1 2 
H at g= 
Bo (40) ee 7 18 


2 1 a ee | 
Signe ee a eee | Se ae 


B, (40) 
i = 4,6,8,...16 
1 (-H,. + HA) -eittk at a=l 
7K}, 1Se 16 20 ee 18 
where 4o = 3 and Bir Hy and Ay are values at i Gaia ie i even. 


Substituting Haq = 2Hyg - Hyg into the last equation one has the 


tri-diagonal coefficient matrix 


ane “eee . _ —s ee See 
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To solve for G, the centered finite-difference scheme is 


applied to the expanded form of Eq. (11); 


2 
1;9°G r) G7 _1l 
7 BLS a 2 “ee a) - co? : 


at interior levels. The discretized equations are 


G 
Shes eG" Ge} e259 atce= te 
Scat. <i 18 c 
2 
6B. 
ly; 4 se ee | eae 
Sigel UG pg 2G Giga grit Sieg) = Ge ate 
i 2 c 18 
(Ac) 
ic 1 a Ps 
~ golw_wo!S15-S17)- aoe Sie! pe 3 ne eee 


where SB; /do is the finite difference form for F5tnB; The first 
and last of the above three equations were established by 
evaluating Eq. (11) at o=l and o=17 respectively, and 
incorporating the appropriate boundary conditions at 

o=0 and o=l (i=18), the latter given by Eq. (13). For the last 


equation, Gig = 2617 a Gig and Gi¢ = (1/2) (G),7+G,5) must be 


used. One has the tri-diagonal coefficient matrix 


Son 


(Aa)* 


4+5065B. 
i 


1 


1 
B) 
~38 
~B. 
L 
a + So 
160 @Aj gl 


3354 7p 1S. 
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Appendix C. Solution of the horizontal modes 
The horizontal modes for given equivalent depth (denoted c) 
and for specified longitudinal wavenumber (m) are determined from 


Eq. (17) which is given as 


(A-vI)K = 0, 


0 f mc 
acos > u 
A= f 0 $ ty and xX =| v 
mc -c_ 9 ()cos¢ y 

acos$ acoso 3¢ 0 


where we suppress the subscript m on the vector elements of x . 
After discretization, the veccor x of unknowns consists of 111 
elements for a 5° latitude increment including polar values. 
Each element in A forms a Square block of size 37. Moreover, 


non~derivative blocks are diagonal. Using the 4th-order finite 


differencing scheme given as 
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a band coefficient matrix with width 5 (zero diagonal except the 
two corners) is formed for a derivative block. More 
specificially, and to see how the boundary conditions are applied 
to the numerical scheme, let us write down the basic finite 
difference formula at any latitude other than the pole; e.g., 


85°N: 


dy) a5~ ~Yl¥95~8¥99*8¥gq-¥75) 

where Y = oy av =5° and subscripts indicate degrees of 
latitude, northern hemisphere. For those grid point values at 
latitudes greater than 90°, one must use the following expression 
which is always true and is thus unlike the given boundary 
conditions which only hold at the poles. If A and y represent 


longitude and latitude, respectively, one can show that 


v(Ayy) = -(-1)"v(aAtn, y) 


Ory) (-1)™y(Atn, y) 


by referring to the Fourier expansion given in Eq. (16) and using 


eim(A+m) = (-l yM_imr 


A negative sign is needed in the v-formula because the reference 


frame changes direction across the pole. For example, 


" 
b 
nn 
t 
i 
=) 
< 
foe) 
Ww 


Yo5 


- a m 
Yos = ( 1) Yes, 
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Owing to the differences of boundary conditions, there are three 


mameguenaneey 


cases to consider. Discussion is focused on the northpole. The 


Situation at the southpole is similar. 


Case lL. m2: all three unknowns vanish at the pole. 


ee 


Finite differencing is needed only between 85°N and 85°S. 


Beginning from the northernmost point the first three grid points 


yield 
BY) = ¥(=Voe +8¥ “8 ,t¥a-) 
sy’ 85 95. "90 80°75 
= ¥( t¥g5-8¥gqth75) +/-: if m odd/even 
sh) ag = 1E “V9 948 -8¥ 46440) 
dy’ 80 90°°%85 75° *70 
= BV 35 ~B¥75*¥79) 
oY). = 1 -Vg5t8Vgq 879 +¥e5) 
ay ’75 85°" "80 70°%65’° 
Hence, the upper left corner of the < = block has the form 
ie: ae Yes 
8 0 -8 1 Yeo 
cY .— (C-1) 
¥75 | * —a Wve 


IW) 


the band coefficient matrix alluded to earlier. For the other 


derivative term, one has, noting that y=a9, 


oe re oy ee 
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+ AV gghOS Pye t8V gg 


COS$ gp COS$.4, 


= y (4V,,-78v chy ————.) 
85 80 COS dye 75 COS Oe, 


pee ee ea. = , 
waseag (VCOS9) | 99 COEF, 5 90°90 EV gsCOS tgs “PV 7gCO8H 75 V9qS95% 79) 


J vey Sees 2 BNge oS ge ana 
85 COS $9 COS ¢ 96 70 cCOSdgg 


Except for ratios of cosines, these equations give the same 


coefficient matrix as W. 


Case 2. m=0: y does not vanish at the pole. 


For the first derivative term one has, 


oy s ms x ' 

By)as ~ Y'-¥95*8¥q9 BY a0 +" 75) | 
j 

ay = = = 


It is seen that there will be an extra column for coefficients of 
V9 added to W of Eq. (C-1), expanding its dimension to , 


35x37. For the other derivative term, finite differencing at the 


pole is needed. To achieve this, we write 
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1 3 = I(vcoss)_ 3 
aosy 3p (VOOSe) = ScasRSy > ok, 


treating z=sing as the independent variaple and p=vcos@ as the 


dependent one. Using a second order Taylor's expansion 


(2 2 
; ap “p 857290 
Pao * Pootrz) 90! 4e57 299)* 7) 90 = 


2 
2 (2,,72,,) 
7 aP _ a“P 80 “90 

go ~ P90*s2)90'%g0"790)*-— 2) 90 ——z-—— 


Seer o : ap, 
and eliminating 322) 90" one can solve for 5Z) 90: 


l-singgy COS$ ge l- “Singge COSbgp 
) ==: ( sing 8 - (——_— -—_—_—_ )v 
3Z 90 Singge-Sindsy T-singss 85 Sings o-Singss I=sindsy Vgo° 


Hence, there will he an extra row for coefficients of vgs and vgq 
added to the corresponding block of case 1, expanding its 


dimension to 37x35. The characteristics of the completed 


coefficient matrix is sketched below. 


ee SAPO AERIN Sr i Rt ALA 


diagonal 


row J07 


Case 3. me=l: only ¥ vanishes at the pole. 


One needs a finite-differenee form for sf at the pole. We 


choose a one-sided, 4th-order scheme 


av, 4 ¥907¥%as _ 1 ¥907%RO 
4690 3 7 KS 3 35 


= V¥(-16¥56*2¥g0) ° 


The other derivative term is first exvanded before bend 


aporoximated by finite differencing: 


' ? rayv } 
aLY-T) 33 (vcos@) i gc = ‘Sry = vtane:¢- 


" 


V("Vge t8V gg 8¥ gg t¥ 76) “Veo tanbac, 


= BYVgg ~(yttandor Wye eWogtyv—- 


The remaining discretized cauations are identical to those in 


6? 


In this’ cafe, there bs a (One seto CoV itycient: of ugg tror 
the Coriolis term and one gust use Che toundary candi tian ueev at 
the pole to shaft this coefiserent to vgg. Wus the polar values 
OO ub Merevt ret bee feet greed. "Ws COFFESPON ED COLLEEN and FUWwS are 
therefore remove-+. ‘The qomplete structure of this coefficient 


Matrix as osiven below, 


Por the finite-elerne it method, cach dependent variable is 


excaidded in oa set of linear space functions, (6), 


oe 


wp) % ae 

upd y = = v ftp! re 
ere De } ; 

boo v 


whe re Byte) are BASiCO fuMmEtionas fur sudsnes of Fat Mepeeee 


defined as 


0 Ost, 
9-5,_ 
tt Oya Osh, 
pF yad : 
e,tared 
} 
o-3 
ne Gk b. sass 
Byeiey eee 
0 Pe ha 


Note that we aqgain drop the = subscrint on elomentsa of the 
vector ob. Substituting Fq. (4.2) into fq. (17) and multiplyina 


by acoss, the system yields 


v 


~avcost.s,u,8, * ausin2etiv. Rk, «4 mecca. = 9 
a 4 - ‘ rg 


» 
c he oe * ss -~ 
avcoeety,G. + are Civ. sql By cose) = Q . 
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Pog the Wane when all boumdigy Values are sero, the discretized 


Sfateit Benes 


, 8 0 3 v« ft a 
a a ee ede de oo «4 v de G (C.3) 
0 8 8 7 T # ou Z 


where a - Cage Ug ape cece dye coos UZe vp, and similarly 

foc voamt >. The Bloens of the coefficient matrix are to be 
Conde ructed as follows. OM multiplics cach term by B.(¢) ani 
wateseates Over the domain. For instance, from the first term 


One Dw ea, 


e*Q % Kel 
: -aveoseR ve) fou Bi(a)de # = 3 $,.uU 
at/e Pe hd Sere yBkel "kj 3 
where 
o oe R4 
2.) 8 apne 46 


.@ the element of S at the kt row and anh column. Notice that 


HR, 20 if -k-zIX]L. Sitilarly, 


-"@e i tee 
f = ; ausinZsR R_ do 
©: ar Roy 
<~) 
a. 
A * 
t. = aes mchR Rod 
Lan a eee < } 
mms 
‘ eel ante 
=O Ro. ee! 
we o, ccosaR, a Ve. 


x-) 


It must be pointed out that 


= 0/2 5 
W,. = -¢ (B.cos¢)d¢ 
kj 1/2 Be wet 
n/2 es ob, 
= -cB, B.cos¢ | + cR.cos¢ ds 
ae. -n/2 -0/2 J 9 
$ aB 
k+l] k > 
Mics ccos¢B. xp 9 = Woe 


This shows that matrices W and W are the transpose of one 
another. Once again, all blocks are square and tridiagonal with 
dimensions of 35 provided that all variables vanish at the 

pole. Furthermore, Eq. (C.3) can be put in the following 


standard form: 


(A-v1)x = 0 (C-4) 
0 g7lp slp 

a ={s71 ¢ 0 s7lw 
sly s7lh 0 


and s7! is the inverse matrix of S. 

For the case of m=0, y does not vanish at the poles. The 
bottom-right S-block of (C.3) is consequently expanded to 37x37, 
W to 35x37 and W to 37x35. One still obtains a standard 


eigenproblem. In the case of m=l, u and v have non-zero boundary 


values. The A matrix is thereby expanded to (374+37435)2, and the 


Sie Oe ai ae - oo 
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details of (C.4) look as fo.lows: 
Iy9x37 9 Q "37 
vf o 139%37 4 Vag * 
0 0 135%35 35 
0 sz F sz T u 
37x37 37x37 37x37 737x35 37 
a : ay 
339437 ° 37x37 0 339x%37 37x35 V37 
ay et + 
5 35x35 735%37 $35%35 35x37 0 Yas). 
Finally, using the boundary conditions of u) = -vy and 37 = V37, 


one can shift column 1 and 37 to 38 and 74, respectively, then 
remove columns and rows of 1 and 37. The resulting system has 


dimensions of 107. 
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Appendix D. The Shooting Method 


Rewrite the differential equations (18) in matrix form as 
y V=EV 


and approximate by an implicit finite-difference scheme such that 


Vv. 


- Vv. Vv. + ¥. 
j+l lees 3+1 j 
BY tT a ice 
or 
= 2 _ 4y ~-1 Ay <= 
= EU¥j41/2i¥)¥5 
where vs = [uly ), viy3)17 and F is a 2x2 matrix with elements 


functionally dependent on Y341/2 and v . In addition, the rate 


of change of v5 with respect to v at each grid point is computed 


(D.2). 


At the end of the integration (y=yj), Vy can be forced to 
converge to Vp (the value specified by the boundary condition) by 


a modified frequency, Vn according to Newton's method. Using a 


numerical approach, 


Pema agi eee 
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ernie (D~3). 
m 


Vig may be calculated from (D-3) and the derivative may be 


calculated from (D-2). Since we only need the boundary value 
(point J), (D-2) is integrated along with . noting only that 
34) 520 vanishes because of the fixed boundary (initial) value. 


To cycle the shooting method, the new value ie) is used in (D- 


1). 
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Figure Captions 


Vertical structures of the vertical motion field. There 
are 9 eigenvectors solved from the H-equation (Eq. 10), 
ordered by decreasing equivalent depth (left to right, 
top to bottom). The top (ordinate level 1) is 

at o = 1 and the bottom o = l. All unmarked abscissas 


range between £1. 


Vertical structures of the wind and height fields. 
There are 9 eigenvectors solved from the G-equation (Eq. 
11), ordered by decreasing equivalent depth. level 1 is 


at o = i level 9 at ae , and interior levels are 
2jtl 

at 18 

Comparison of modes for 5° and 10° latitudinal 


increment. 
Boundary modes for channel model between 60°N and 60°S. 
Low-order Rossby modes as they depend on model geometry. 


Eigenfrequencies for longitudinal wavenumbers m = 0 and 
1 with equivalent depth 9.555 km (external modes). The 
abscissa is the order of magnitude of frequency in 


sec7!. the ordinate is the order of frequencies from 


the most negative (#1) to the most positive (#107). 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 
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First and last 35 frequencies represent westwerd and 
eastward gravity modes, respectively. The last mode 
with negative frequency is #58. There are no values in 
the range of Rossby modes from #36 to #72 for m = 0. 
Elsewhere, if neighboring values are too clos«, only one 


of them is plotted. 


Same as Fig. 6 except for m = 2, 6, 10. Zero's are used 
for m = 10. 
Same as Fig. 6 except for m=6. Frequencies for all 9 


equivalent depths are plotted (finite difference). 


Same as Fig. 8 except that the finite element method is 
used to establish the eigenproblem. Notice that there 
are more negative and relatively less positive Rossby 


modes than for the finite difference result. 


Latitudinal eigenstructures of u for m = 6 and 
equivalent depth 9.555 km. Shown are westward gravity 
modes from #21 to #35, each plotted from pole-to-pole at 
a 5° increment. Whether a mode is computational or 
Physical is denoted by ac or a p and they are 


sequenced. 


Same as Fig. 10 except for results from the finite 
element method, 


Amplitudes of real data set projections onto the 


Fig. 


Fig. 


13 


14 
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latitudinal modes of the 4-th order difference model for 
m = 6 and equivalent depth 9.555 km. The abscissa 
represents order of magnitude of amplitude and the 
ordinate gives order number of modes. Physical modes 


are indicated by connecting segments. 


Same as Fig. 12 except for the finite-elemenmi odel 


modes. 


Comparison of original and reconstructed latitudinal 
profiles. Dotted line is the real part of observed u 
after projecting onto the first (external) vertical mode 
for m=6. Solid line is the reconstructed profile by 
removing amplitudes of computational modes using the 
finite difference method for normal mode analysis. 


Broken line is the reconstructed profile based on 


filtering with the finite element modes. 
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Fig. 1 Vertical structures of the vertical motion field. There are 9 ciaen~ 
vectors solved from the H~cquation (kqg.10), ordered by decreasing 
} equivalent depth (left to right, top to hottom). The top(lordinate 


level 1) is at a = z and the bottom 9 = 1. All unmarked abscissas 


range between ‘1. 
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Fidq. Vertical structures of the wind and height fields. There are 9 


to 


eigenvectors solved from the G-equation (ka. 11), ordered by 
: : 7 
decreasing equivalent deptn. Level 1 is at a = +5 level 9 at a 
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and interior levels are at are 
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BOUNDARY MODES 
(CURVATURE , 60°N - 60°S ) 


fig. 4 Boundary modes for channel model between 60°N and 60°S. 
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Has Vaca pienfre quencies for lenitudinal wavenumbers m= % and 1 with cquivalent 
fh of: tao km fexternal mode ). The abscissa is the order of macnitudce 


-4 : 
Totremence in see oo. Th ordinate is the order of frequencies from the 


mot negative (#1) ta tte. wont oo ative (@107).  Firet and last 35 
Sreqencte > repres:ut westeard and «eastward gravity modes, respectively. 
The Tact modes with neqative fire caency is &'&. There are ne values in 


the range of kowshy modes from #36 to #72 from = 0. Elsewhere, if 


noighboring values are tor clos, only one of them is plotted. 
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fie. & Same as Fig. 6 exce;t for m = 6. Frequencies for all % equivalent depths 


are plotted (finite differsnce). 
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Fig. 9 Same as Fig. 8 except that the finite element method is used to establish 
the cigenprohlem. Notice that there are more neaative and relatively 


less positive Rossby modes than for the finite difference result. 
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Amplitudes of real data sect projections onto the latitudinal modes of 
the 4th order difference model for m = 6 and equivalent depth 9.555 km, 
The abscissa represents order of maqdnitude of amplitude and the ordinate 


gives order number of modes. Physical modes are indicated by connect- 


ing segments. 
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finite-element model modes. 
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Fig. 13 Same as 
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Compar isin of igina) and Yeconstructog latitudina] Profiles, Dotted 
line js the roa Part of observed u after Projecting Onto the first 
(oxterna}) vertical mode for m = 6. Solia line is the reconstructeg 
PYofile by removing aM li tudos of computationa) Modes Using tho 

finite di f€oronee Method for normal mody analysis, Broken line is 

the Feconstructed Profile based on filtering with the finite eloment 
Modes, 


